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ABSTRACT 

We study the eccentricity distribution of a thick disc sample of stars observed in 
the Radial Velocity Experiment (RAVE) and compare it to that expected in four 
, simulations of thick disc formation in the literature (accretion of satellites, heating 

$_j ■ of a primord ial thin disc durin g a merger, radial migration, and gas-rich mergers), as 

compiled by Sale s et ail (|2009| ). We find that the distribution of our sample is peaked at 
low eccentricities and falls off smoothly and rather steeply to high eccentricities. This 
distribution is fairly robust to changes in distances, thin disc contamination, and the 
particular thick disc sample used. Our results are inconsistent with what is expected 
for the pure accretion simulation, since we find that the dynamics of local thick disc 
stars implies that the majority must have formed in situ. Of the remaining models 
explored, the eccentricity distribution of our stars appears to be most consistent with 
the gas-rich merger case. 

Key words: Galaxy: disc - solar neighbourhood - Galaxy: formation - Galaxy: 
structure 



1 INTRODUCTION 

The thick disc has been a known component of the 
Milky Way for over 20 years (|Yoshiilll982llGilmore fc Reidl 
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1983). Analogous components have been identified in ex- 
ternal galaxies, re vealing that th i ck discs may be quite 
generic features dBursteinl 1 19791 : Ivan der Kruit fc Searld 
Il98ll : lYoachim fc Dalcantonl (20061 ) . Most of the thick disc 
formation scenarios that have been proposed fall into one 
of the following four categories: accretion, heating via a mi- 
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nor merger, radial migration, and intense star formation in 
a gas-rich turbulent environment. 

In the accretion model, satellites infall on coplanar or- 
bits and form the thick disc as they are disrupted and in- 
corporated into the main galaxy; in this case, the thick disc 
would thus consist of stars ori ginating largely (m ore than 70 
per cent in the simulations of lAbadi et al.| [2003) in several 
disrupted satellites. 

The most often discussed scenario for the formation of 
a thick disc is the accretion of a massive satellite by a pre- 
existing disc galaxy, which is thus heated dynamically. The 
resulting thick disc is mainly made from stars that origi- 
nated in the primary galaxy's primordial thin disc rather 
th an the merging sate l lite. F or example in the simulations 
bv lVillalobos fc Helmil (|2009h . 10 to 20 per cent of the stars 
in the remnant thick disc at the "Solar" radius are accreted. 

During a turbulent gas-rich phase, stars may also form 
in a thick disc. This can happen in massive clumps in a 
rotationally supported component, as in the simulations of 
iBournaud et al.l d200E ), but also during gas-rich mergers as 
shown by Brook et al (|2004h . In the latter case, the result- 
ing thick disc stars would have been formed in situ and with 
relatively hot kinematics. 

Merging events of any variety might be unnecessary to 
explain the phenomena, however. Stars may migrate radially 
from the inner parts of a galaxy to the outer regions due to 
reson ant interacti ons with spiral arms ([Sellwood fc Binnevl 
I2OO2I ) and a bar l|Minchev fc Famaevl 120091 ) . Although the 
migration process itself does not heat the disc, a greater 
vertical velocity in the higher surface brightness central 
regions results in larger heights above the plane being 
reache d in the outer regions where the surface brightness is 
lower (|Kregel. van der Kruit. fc Freeman] 120051 ) . Migration 
can thus result in the formation of a t hick disc from thin 
disc stars without any ext ernal stimulus (|Roskar et al.ll2008l ; 
ISchonrich fc Binnevl [20091 ) . 

Recently. [Sales et alj l|200gh suggested a dynamical test 
that could differentiate between the formation models and 
be applied to data to disentangle which is the most likely to 
have occurred in our galaxy: they propose using the eccen- 
tricities of thick disc stars as a discriminant, as populations 
formed in situ are likely to move on low eccentricity orbits, 
while those accreted can have any eccentricity, but will typi- 
cally be biased towards higher values. As kinematic data for 
thick disc stars become available fr om surveys such as th e 
Radial V elocity Experiment (RAVE Isteinmetz et afll2006l) . 
SEGUE l|Yannv et all 120091 1. and eventually Gaia, we can 
apply these tests to see what they can reveal about thick 
disc formation. 

In this paper, we investigate and constrain the dynamics 
of a sample of Milky Way thick disc stars and compare the 
resulting eccentricity distribution to what is expected for the 
above four models. We use data from the RAVE survey and 
distances calculated in the manner of iBreddels et al l l|'201Ch 
but with modifications described in Section [2] In Section [3] 
we derive the eccentricity distribution of the stars and in Sec- 
tion|?]we then examine this in light of the accreti on, heating, 
merge r, and migration simulations discussed in lSales et al.l 
|2009l ). 



2 DATA 

2.1 The RAVE survey 

ravbH measures radial velocities and stellar atmospheric 
parameters from spectra using the 6dF multi-object spec- 
trometer on the Anglo- Australian Observatory's 1.2 m UK 
Schmidt Telescope. The survey looks in the Ca-triplet region 
(8410-8795 A), has a resolution of ~ 7500, and is magnitude 
limited. The targets chosen are southern hemisphere stars 
taken from the Tycho-2, SuperCOSMOS and DENIS sur- 
veys with I-band magnitudes between 9 and 13. The average 
internal errors in radial velocity (RV) are ~ 2 km s _1 , and 
the approximate RV offset between RAVE and the litera- 
ture is smaller than ~ 1 km s _1 . The catalogue also includes 
2MASS photometry and proper motions from Starnet 2.0, 
Tycho-2, SuperCOSM OS, and UCAC2. Fo r more informa- 
tion about RAVE, see lZwitter et all (|2008f ). 

2.2 Distances 



Distan ces were calculated along the lines of Brcddcls et al. 
(2010). We briefly sketch the method here and refer the 
reader to IBreddels et all (|201Gh for more details. 

The stars are fit u sing the Y 2 (Yonsei-Yale) isochrones 
l|Demaraue et al]|2004l ) and their measured characteristics. 
The measured quantities of the stars (T eg, log(g), [M/H], J , 
and J - K s ) from the RAVE pipeline (IZwitter et al .1120081 ) 
and 2MASS are used to minimise the x 2 statistic to find 
the closest model star for each observed star. Then, the er- 
rors of the observed quantities, which are assumed to be 
Gaussian, are utilised in a Monte Carlo simulation, from 
which the absolute magnitude and its error are determined 
from the resulting probability distribution function. Stars 
for which none of the isochrones provide adequate matches 
are discarded. Since the Y 2 isochrones do not extend past 
the RGB tip, clump stars may result in poor fits; in addi- 
tion we have explicitly removed the clump stars that should 
remain in the data set (see Section f2 . 3 p - 

The effecti veness of this distanc e method was exam- 
ined in detail in IBreddels et ail (|2010h . but the main points 
will be mentioned here. Comparison to main sequence 
stars observed by HIPPARCOS showed good agreement in 
the measured paralla xes. To test the distances for giants, 
IBreddels et al.1 (|201dh used the members of M67, an old open 
cluster. The mean of the calculat ed distances was 1.48 ±0.3 6 
kpc, compared to 0.78 kpc from Vand enBerg et al. I (|2007l ). 
This discrepancy is related to the fact that when age is left as 
a variable, stars on the giant branch are sometimes fit better 
by (unrealistically) younger isochrones. However, when us- 
ing the 4 Gyr isochrones, which corresponds to the accepted 
age for the cluster, the distances agreed with those in the 
literature within error (see Figure [TJ . 

This has motivated us to calculate distances for all stars 
setting the ag e at 10 Gyr, which is th e characteristic age of 
the thick disc fEdvardsso n et al Il993l ). This implies that our 
method will assign incorrect distances to younger thin disc 
stars, but since younger giant branches are brighter than 
older ones, this assumption results in the assigned distances 

1 We use the 30 August 2008 internal data release, which consists 
of 135,338 stars. 
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Figure 1. The theoretical solar metallicity, 4 Gyr isochrone with 
the CMD of M67 giants. Points with crosses are the M67 stars 
transf ormed to absolute magnitudes using the iBreddels et al.l 
d2010h distance calculated after excluding the clump stars of 1.48 
kpc; filled points use the distance of 0.78 kpc calculated assuming 
an age of 4 Gyr. Gray points are stars identified as belonging to 
the red clump. 



to younger stars being slightly smaller than they actually 
are. As discussed in Sec. 12.31 we will be selecting stars with 
1 kpc ^ \z\ ^3 kpc, so as to isolate a thick disc sample. The 
consequence of using the "wrong" isochrone for young stars 
is to reduce their \z\ and to move many from lower \z\ out 
of our sample rather than scattering thin disc stars up into 
it. Even so, we investigated what effect that would have on 
the distances and the level of thin disc contamination more 
quantitatively in Section f3. 2. 21 



Figure 2. Velocity distributions of our thick disc sample (1273 
stars) computed in a right-handed reference frame at rest with 
respect to the Galactic centre. Mean velocity values and disper- 
sions, calculated using robust estimation, for each component are 
given in each panel. Stars to the left of the dashed line in V y were 
discarded when calculating the means and dispersions. 

RAVE's metallicity pipelin fl but also because this mimics 
more closely the selection bv lSales et al.l (|2009l ) . We further 
clipped our sample, discarding all stars with V y < 50 km s 
to minimise contamination from the halo, and including 
only stars within a heliocentric cylinder with a radius of 
3 kpc so the data was in a form best suited for comparison 
with the eccentric ity distributions of models as illustrated by 
ISales et al.l l|2009l ). The final sample consisted of 1273 stars. 
The velocity distributions are given in Figure(2] The velocity 
dispersions are generally consistent with the literature val- 
ues, especially in V x , while those in V v and V z are slightly 
highe r by approximately 6 — 8 km s _1 (jSparke fc Gallagher! 
12006! ) . 



2.3 Sample selection 

We first cleaned the data set by discarding any stars with 
distance errors > 40 per cent, proper motion errors in ei- 
ther RA or DEC > 10 mas/yr, or radial velocity errors 
> 5 km s _1 . Once the age was set at 10 Gyr, the clump 
stars were not fit well by any of the isochrones, so most of 
them should have been discarded on that basis. To ensure 
all clump stars were indeed removed, however, we threw out 
stars with log(g) > 1.5 and J — K s < 0.75. 

In order to isolate a sample of thick disc stars, we 
chose the remaining stars with \z\ between 1 and 3 thick 
disc scale-heights, which corresponds to the range 1-3 kpc 
ijVeltz et al] 1200a ). The decision to make our thick disc se- 
lection based only on \z\ rather than including a metallic- 
ity criterion was motivated partly by uncertainties in the 



3 RESULTS 

3.1 Eccentricity distribution 

To calculate the eccentricities of the RAVE stars in our 
thick disc sample we integrate d their orbits in a Galac- 
tic potential. This consisted of a lMivamoto fc Nagail (|l975l ) 
disc, a Hernquist (1990a) bulge and a spherical logarithmic 
halo. In this model, the characteristic parameters used were 
Mdisc = 8.0 x 10 10 ,M bu igo = 2.5 x 10 10 ,^ a io = 27000.0, a = 
6.5, b — 0.26, c = 0.7, and d — 12.0, with masses in Mq, 
velocities in km s _1 , and lengths in kpc, which produce a 

2 For completeness, we have tested that our results do not change 
appreciably when we focus only on the subset of stars with high 
S/N spectra which according to their [M/H] belong to the thick 
disc. 
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circular velocity of ~ 220 km s at 8 kpc from the Galac- 
tic centre. The eccentricities of the stars were defined as 



(r. 



ri)/(r apo + r pcr i), where r apo (r pcr i) is the maxi- 



mum (minimum) distance reached by the star in its orbit. 
Figure [3] shows the eccentricity distribution obtained for our 
sample. The main features of this distribution are an asym- 
metric peak at a fairly low eccentricity value with a rel- 
atively (though not entirely) smooth falloff towards higher 
eccentricities. This peak at low eccentricity suggests that the 
thick disc was formed primarily in situ (see Section [4.2[) . 



3.2 Robustness of the eccentricity distribution 

3.2.1 Systematic errors 

We expect that the largest contribution to systematic er- 
rors would be due to the distances. Distance overestimation 
would result in larger heliocentric velocities, because the ob- 
served proper motions would be placed at larger distances. 
Thus, the distance overestimation should result in a peak 
at a higher eccentricity than more accurate distances would 
reveal, not a peak at lower eccentricities, suggesting that 
the calculated distribution would not likely be a result of 
distance overestimation. On the other hand, distance un- 
derestimation should result in smaller heliocentric velocities 
and more circular orbits, and thus could cause an artificially 
strong peak at low eccentricity. 

To explore more quantitatively how errors in distance 
could affect the final eccentricity distribution, we calculated 
eccentricities for the thick disc sample using distances that 
were larger and smaller by 20 percent. This value was chosen 
because nearly all the stars in the sample had distance errors 
smaller than 20 per cent, with a peak at 10 percent, making 
this figure conservative. As Figure shows, the distribution 
with smaller distances has a slightly higher peak at a slightly 
smaller eccentricity than the original distribution, and that 
with larger distances has a slightly lower peak at a slightly 
larger eccentricity, as expected. However, both distributions 
remain reasonably close to the original one, which suggests 
that a systematic distance error of 20 percent would not 
result in a substantial change in the thick disc's eccentricity 
distribution. We have found that even with up to 40 percent 
larger distances, the peak in eccentricity remains below 0.4 
and the generally triangular shape does not change. 

We have also tested how our proper motion errors af- 
fect the eccentricity distribution. We found that the main 
effect is to slightly lower the amplitude of the peak at ec- 
centricity ~ 0.2 and to marginally increase the number of 
stars with high eccentricity. Therefore we conclude that the 
shape of the eccentricity distribution is generally robust to 
the estimated uncertainties in our observables. 



3.2.2 Investigation of possible thin disc contamination 

Since the thin disc characteristically has stars on fairly cir- 
cular orbits, there is a possibility that the peak at low ec- 
centricity is caused partially by thin disc contamination in 
our sample. 

To give a quantitative estimate of this contamination, 
we formulated a simple model. Using the Padova isochrones 
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Figure 3. Eccentricity distribution of the thick disc sample using 
the original distances (1273 stars, solid histogram) and distances 
that are 20 percent smaller (1350 stars, dashed histogram) and 
larger (1204 stars, dashed-dottcd histogram). The differing num- 
bers of stars are a result of discarding stars with low rotational 
velocity, since this number depends on what distance is assumed. 



and the IChabrierl i|200ll ) initial mass function (IMF), we 
generated two samples of stars. A thin disc population was 
created from a solar metallicity, 5 Gyr isochrone, and a thick 
disc population was created from an isochrone with [Fe/H] 
= —0.6 dex and an age of 10 Gyr. The relative fraction of 
thin to thick disc stars /thm2thick was calculated using the 
expression 

. -M/ z thin+W/ z thick 



,/i 



where / nolm is the relative fraction at the Sun, and Zthin and 
zthick the thin and thick disc scale heights respectively. We 
ass umed Zthm = 225 pc and z t hick = 1048 pc, as estimated 
bv lVeltz etail (120081 ) for the RAVE sam ple. The local nor - 
malisation / norm was calculated using the lVeltz et~all (2008) 
ratio of thin to thick disc dwarfs and the Padova isochrones 
to determine the fraction of dwarfs for the thin and thick 
disc samples. We assigned z from 1050 pc to 2950 pc in in- 
crements of 100 pc and selected the model stars that would 
have been observed by RAVE (9 < / < 13). To synthesise 
the effect of calculating distances assuming an age of 10 Gyr 
on the younger thin disc population, we found the absolute 
magnitude on the giant branch that would have been as- 
signed to each thin disc star using a solar metallicity, 10 
Gyr isochrone based on their effective temperatures. New 
distances were calculated for the stars with the new abso- 
lute magnitudes and the original apparent magnitudes and 
then re-binned based on these new distances. 

Ratios of the number of observed thin disc stars (both 
using the 5 and 10 Gyr ages) to observed thick disc stars 
were then computed for each bin in distance and overall. In 
the lowest z bin, the thin disc contamination was as large 
as 30 percent, but it dropped to about 10 percent for the 10 



Testing formation mechanisms of the Milky Way's thick disc 5 



7-71 Original Sample 

- 1 / > 



XX X 
X\X 



1 kpc < |z| < 1.3 kpc 
1.3 kpc < |z| < 3 kpc 



v y >. v \ 

< XX >< K7< 

. X A A, A X 

: ;< a x X \ X 

* X > X A 
X > X \f >■ v 



A X y\ X ,'< A J\ -\- \ 

X A x X x X X H, 
'/ v V > V V V V ,\ 

^XXX/vKXA > 

X / X/ N v N ^ x ^ N ./ N J X \ 

/ \ ✓ \ ✓ \ ✓ ^ ^ / h \ - 

A .A.. /\ /\ /\ /\. /\ /\ 7\ y\ 



X X X v y v y y y y y ■ 
N \ V\ ^\ ^\ >\ ^\ ^\ 'X ^\ ^ 

^ A A A A A A A A A A A A 
A A A >f A A A 'x A. A A 'A. x / y y y y _ ( >- l 

O 0.2 0.4 0.6 0.8 

Eccentricity 



^0.2 
o 

^ 0.1 



i^O.2 
o 

^ 0.1 



1 1 1 1 1 

~<V x > = -28.8 +/- 0.5 km/s 
" <T X = 64.3 +/- 0.4 km/s 

■<V X > = -16.8 +/- 0.5 km/s " 
" (7 X = 58.6 +/- 0.5 km/s 


<V > = 166.9 +/- 0.5 km/s " 
" CT^=47.4 +/- 0.4 km/s 

■<V'> = 176.1 +/- 0.6 km/s " 
' CT =43.2 +/- 0.4 km/s 


; v A : 




1 ^Trf 1 1 1 | 1 1 f~j-i-l 1 | 
"<V Z > = -1.8 +/- 0.5 km/s ■ 
" CT Z = 50.0 +/- 0.4 km/s 

"<V Z > = -4.1 +/- 0.5 km/s 


-400 -200 200 400 
Velocity (km/ s) 


" CT z = 39.0+^^^^ta/s 





-400 -200 200 400 
Velocity (km/s) 



Figure 4. Eccentricity distributions of the original sample (solid), 
the subsample of stars with 1 kpc < \z\ < 1.3 kpc (dashed his- 
togram, 450 stars), and the subsample with 1.3 kpc < \z\ < 3 kpc 
(dashed-dotted histogram, 823 stars). 

Gyr thin disc by z = 1250 pc and by z = 1350 pc for the 
5 Gyr thin disc. Since the effect of calculating the distances 
of younger stars with 10 Gyr isochrones is to underestimate 
their distances, the thin disc contamination was lower at 
most z intervals for the 10 Gyr thin disc ratios, since moving 
the stars down in z took many out of the z range considered. 
The overall thin disc contamination for the 5 and 10 Gyr thin 
discs were, respectively, 4.7 and 3.1 percent. 



Figure 5. Vel ocity distributions of the thick disc sample using 
the distances of lZwitter et alj l2010h (black) and the clump thick 
disc sample (gray). The mean velocity values and dispersions 
(calculated using robust estimation) after removing stars with 
V y < 50 km s _1 are given in each panel. The Gaussians plotted 
have the same mean and dispersion as the data. 

amount of higher eccentricity stars was noticeable. These 
trends are natural consequences of increased contamination 
by halo stars. Overall, the higher \z\ distribution retained 
the general properties of the original distribution. Thus we 
conclude that the thin disc contamination is not likely to 
have a large effect on the overall shape of the eccentricity 
distribution. 



Since those ratios suggest that there could be significant 
contamination in the portion of our sample that is closest 
to the plane of the disc, we took our thick disc sample and 
divided it into a low \z\ and a higher \z\ portions to see how 
robust the shape of the eccentricity distribution was to thin 
disc contamination. Two separate trials were performed. In 
the first, the division was placed at \z\ = 1.3 kpc, because the 
simple model suggests that the contamination levels have 
dropped below 10 percent by that height when the stars all 
were assigned 10 Gyr ages, as is the case for the data set. 
The resulting eccentricity distributions for each subsection 
were then calculated and are plotted in Figure [4] The dis- 
tribution of the lower \z\ portion is strongly peaked at low 
eccentricity with fewer stars at higher eccentricities, which 
would be expected for a sample contaminated by thin disc 
stars, which have predominantly circular orbits. The higher 
\z\ sample retained the roughly triangular shape of the orig- 
inal sample, however. The peak did shift slightly to higher 
eccentricity, but not significantly. 

We also performed another trial and cut the sample at 
\z\ — 1.5 kpc. The resulting eccentricity distributions of the 
two subsamples were similar to those of the first trial. The 
distribution of stars in the higher \z \ subsample was lumpier, 
since there were fewer stars in it, and a slightly increased 



3.2.3 Comparison of the original eccentricity distribution 
with those using different thick disc samples 

Since the distances are a key component in calculating 
the eccent ricity distribu t ion, w e also used another set of 
distances. IZwitter et al.l l|2010l ) calculated distances for a 
RAVE data set (based on a later internal data release con- 
sisting of 332,747 stars) by means of a photom etric parallax 
method based o n that of Breddels et alj (|2010l ) but using the 
Y 2 , Dartmouth (|Dotter et al.ll200&f ). and Padova isochrones 
uniformly, rather than logarithmically spaced in time and, 
when picking a best match model star, weighting the model 
star picked based on mass. No age constraint to tailor the 
distances to a thick disc sample was applied in this set of 
distances. Using the same cleaning and thick disc selection 
criteria as before, we selected a sample from this catalogue 
(that based on the Y 2 isochrones) and calculated its eccen- 
tricity distribution. This thick disc sample consists of 6173 
stars. See Figure for the velocity distributions. 



In order to have a more independent check on the ec- 
centricity distribution, since the above two methods of cal- 
culating distances are quite similar, we selected a thick disc 
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Figure 6. Eccentricity distributions of thick disc samples using 
a) the original distances, b) the Zwitter et al. distances, and c) 
the clump sample. 
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Figure 7. Comparison of the eccentricity distributions of each 
thick disc formation model for stars in the range 1—3 (thick disc) 
scale-heights and cylindrical distance 2 < R/Rd < 3. 



sample of clump stars from RAVE. To choose this sample, 
we used the same cleaning criteria on the full data release 
but did not impose the restrictions aimed at eliminating the 
clump stars. We identified the clump by colour and gravity, 
taking it to have 0.6 < J - K s < 0.7 and 1.5 < log(g) < 
2.7. Then, we calculated the distances using the apparent 
magnitudes and ass uming an absolute c l ump magnitude of 
Mk = —1.61, as in lAlves fc Saraiedinil (| 19991 ). Finally, we 
imposed the same restrictions in |z|, volume, and velocity 
as before and calculated the eccentricity distribution. The 
resulting sample includes 3573 stars, and the velocity distri- 
butions are given in Figure [S] 

The velocity distributions i n Figure [5] show good a gree- 
ment with those based on the iBreddels et all (2010) dis- 
tances. As expected, the velocity dispersions for the red 
clump sample are slightly lower, since this sample should be 
devoid of halo contamination (as red clump stars are only 
present in young or intermediate age populations). Figure [S] 
shows that also the thick disc eccentricity distributions agree 
well with each other. All are strongly peaked at low eccen- 
tricities, have a generally triangular shape, and do not have 
a secondary peak at hi gh eccentricity. The sample based 
on IBreddels et all l|2010l ) falls off at high eccentricities less 
smoothly than the other two, but that could be due to sta- 
tistical fluctuations in this much smaller data set. On the 
other hand, the clump sample has its peak at a lower eccen- 
tricity and contains a smaller number of stars at higher ec- 
centricities as expected since this region is populated mainly 
by halo stars which are only present in the red giant sam- 
ples. So while the three distributions are not identical, they 
agree fairly closely. Results of the application of KS tests 
show that all three are consistent with being drawn from 
the same distribution. This leads us to conclude that the 
eccentricity distribution is reasonably robust to the exact 
thick disc sample selection. 



4 COMPARISON WITH THEORETICAL 
MODELS 

4.1 Discussion of the models 

We now compare our calculated eccentr icity distribution 
to those computed by ISales et al.l (|200"9l) from th e simu - 
lations of lAbadi et all (I2003TI. IVillalobos fc HelmH l|2008h . 
IRoskar et all l|200ct ). and iBrook et al l l|2004 ). These simu- 
lations have been discussed in the literature, so we will only 
briefly describe them here. For d etails, we re f er the reader 
to the literature and to Table 1 of ISales et al.l l|2009l) . which 
sum marises their main parameters. 

lAbadi et al.l (|2003h demonstrate the formation of the 
thick disc through accretion of satellites during the hi- 
erarchical formation of a Milky Way-like galaxy in the 
A cold dark matter paradigm and using cosmological TV- 
body /smoothed particle hydrodynamics (SPH) simulations. 

In the h e ating scenario simulation of 
IVillalobos fc Helmil l|2009l ). a satellite merges with a 
primary thin disc (of comparable mass) on a prograde orbit 
inclined 30 degrees. The subsequent formation of a new thin 
disc from cooling gas and any changes that might cause in 
the thick disc were not modelled in t his sim ulation. 

The simulation of IRoskar et al.l l|200Sl ) models the for- 
mation of a galactic disc by starting with a dark matter halo 
and a hot gas halo which over 10 Gyrs cools and forms stars 
in a disc which migrate as transient spiral structure forms. 

IBrook et alj (|2004l ) model the formation of a disc in a 
semi-cosmological TV-body/SPH simulation which incorpo- 
rates gas heating and cooling, star formation, feedback, and 
chemical enrichment. Their thick disc develops during gas 
rich mergers. 
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ISales et all l|2009l ) found that the location of the peak 
and the shape of the eccentricity distributions are driven by 
what kind of stars make up the thick disc in each model, 
as shown in Figure [7] In the accretion scenario, the thick 
disc is composed mostly of accreted material; since satel- 
lites tend to be on radial orbits, the resulting eccentricity 
distribution is relatively broad and peaks in the mid to high 
eccentricities. The other three scenarios involve a thick disc 
being composed primarily from stars formed in situ. Such 
stars have more circular orbits, which produce a peak at 
low values of the eccentricity. The heating and the gas-rich 
merger simulations have a contribution of accreted stars to 
the thick disc, which show up in the eccentricity distribution 
as a distinct second but much smaller peak at high eccen- 
tricity in the heating scenario and a lumpiness at mid to 
high eccentricities in the merger case. 

4.2 Comparison with the models 

Comparison of Figures [6] and [7] shows that the eccentricity 
distribution of our sample does not support the accretion 
model. The triangular shape with a peak at low eccentrici- 
ties of our distribution does not resemble the broad mound 
peaked at middling eccentricity that characterises the ac- 
cretion scenario. The eccentricity distribution of our thick 
disc sample with its prominent peak at low eccentricities 
is consistent with the distribution displayed by the heat- 
ing, migration, and merger scenarios. Our distribution does 
not exhibit a secondary peak at high eccentricities like that 
evident in the heating model due to the accreted satellite. 
However, the location of this secondary peak depends on the 
initial orbital configuration of the accreted satellite. If the 
initial conditions were changed to produce a secondary peak 
at middling eccentricities this could probably be obscured in 
the data by the larger fraction of in situ stars. 

While the distributions from the migration and merger 
models are not identical to the distribution of our sample, 
they do show some resemblance and have no features that 
seem to count against either being possible. The gas-rich 
merger's distribution even displays an asymmetric peak like 
our distribution, making this scenario the most consistent 
with our data. The migration distribution exhibits a sym- 
metry about the peak (is more Gaussian-like) until it gets 
down into the high eccentricity tail that is not displayed in 
our distribution, making the migration scenario somewhat 
less consistent with our data. 

There are several reasons for the models not being a 
perfect match to the data. Random and systematic errors in 
the data, while not significantly affecting the overall charac- 
teristics of the distribution, can alter the exact shape, espe- 
cially at high eccentricities. Also, since the simulations did 
not attempt to duplicate the Milky Way exactly, they might 
not be similar enough to the Galaxy to precisely mimic the 
eccentricity distribution even if the formation mechanism is 
correct. It is also plausible that the thick disc was formed 
by a combination of processes. Radial migration has been 
shown to be dynamically possible in a galaxy such as our 
own; mergers, including gas-rich, are known to occur. Thus, 
the nature of the measured eccentricity distribution could be 
indicating that both radial migration and gas-rich mergers 
contributed to its formation. 

One possible way to investigate which model is the most 



likely, would be to compute the eccentricity distributions of 
the models and the data for different locations along the 
Galactic disc, as one may expect the contribution of accreted 
populations to change with distance and become more domi- 
nant in the Galaxy's outskirts. Undoubtedly, data from Gaia 
when it becomes available will aid in clarifying how exactly 
the Galactic thick disc formed. 



5 SUMMARY AND CONCLUSIONS 

We have isolated a sample of thick disc stars from RAVE 
survey data, calculated its eccentricity distribution, and de- 
termined that the distribution is fairly robust to changes 
in distances, thin disc contamination, and the specific thick 
disc sample used. Our eccentricity distribution is fairly trian- 
gular in shape, depicting a dominant peak at low eccentricity 
and a relatively smooth falloff at high values. 

We have c ompared th i s find ing with the eccentricity 
distributions in ISales et al. (2009) presented for simulated 
thick discs formed via accretion, heating via a minor merger, 
radial migration, and gas-rich mergers. The broad peak at 
moderately high eccentricities of the accretion model is not 
consistent with the relatively narrow peak at low eccentric- 
ity displayed by our sample. This indicates that the Galactic 
thick disc formed predominantly in situ. A lack of a distin- 
guishable secondary peak at high eccentricity further sug- 
gests that if any of the thick disc was accreted, its direct 
contribution of stars in the Solar neighbourhood was min- 
imal. The gas-rich mergers simulation, and to a somewhat 
lesser extent, the radial migration model, are consistent with 
the distribution of our sample. This suggests that these for- 
mation mechanisms could have had some role in the forma- 
tion of the Milky Way's thick disc. 
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